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Abstract. The decoherence of a two-state system coupled with a sub-Ohmic 
bath is investigated theoretically by means of the perturbation approach based 
on a unitary transformation. It is shown that the decoherence depends strongly 
and sensitively on the structure of environment. Nonadiabatic effect is treated 
through the introduction of a function which depends on the boson fre- 
quency and renormalized tunneling. The results £1X6 cLS follows: (1) the non- 
equilibrium correlation function P(t), the dynamical susceptibility x"^) an d 
the equilibrium correlation function C(t) are analytically obtained for s < 1; 
(2) the phase diagram of thermodynamic transition shows the delocalized- 
localized transition point a\ which agrees with exact results and numerical 
data from the Numerical Renormalization Group; (3) the dynamical transi- 
tion point a c between coherent and incoherent phase is explicitly given for the 
first time. A crossover from the coherent oscillation to incoherent relaxation 
appears with increasing coupling (for a > a c , the coherent dynamics disap- 
pear); (4) the Shiba's relation and sum rule are exactly satisfied when a < a c ; 
(5) an underdamping-overdamping transition point a* exists in the function 
S(u>). Consequently, the dynamical phase diagrams in both ohmic and sub- 
Ohmic case are mapped out. For A <C u c , the critical couplings (ai,a c and 
a*) are proportional to A 1_s . 
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I. Introduction 



A two-state system (TSS) coupled to a dissipative environment offers a unique testing 
ground for exploring fundamental physics of quantum mechanics behaviors such as 
tunneling and decoherence [TJ [2] . It also provides a paradigmatic standard model for 
studying quantum computation [3j. An interesting problem in this type of system is 
decoherence resulting from the influence of an environment, which is usually repre- 
sented by an infinite set of harmonic oscillators. This system can be described by the 
spin-boson model (SBM), which reads 

H = ~ Aa x + ]-ta z + u k b\b k + \Y1 9k{b\ + h)a z , (1) 

~ " k k 

Here a x and a z are the usual Pauli spin matrices, A is the bare tunneling matrix 
element between the two states, e describes the bias of the system. Throughout this 
paper we set h — 1 and k B = 1. The quantity g k represents the coupling strength 
of TSS to the kth oscillator. Without dissipative environment, the TSS exhibits 
coherent tunneling between the two states, while the coupling to the environment 
will result in the lose of its phase coherence. 

Impossible as an exact solution of SBM is, both the equilibrium and non-equilibrium 
dynamics in this model have been widely studied by using the Path-integral formal- 
ism, variational method, real time quantum Monte Carlo, numerical renormaliza- 
tion group (NRG), flow equations, etc [1-26]. Much of the physics underlying the 
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model has not been revealed by numerical methods, such as the coherent-incoherent 
transition P, [T7j. The main interest is to understand how the quantum environment 
influences the dynamics of TSS, in particular, how dissipation destroys quantum co- 
herence. On the one hand, when the system can be prepared in one of the two states 
by applying a strong bias for times t < and then let it evolve for t > in zero bias, 
the non-equilibrium correlation function P(t) is of primary interest [U EI] . Moreover, 
P(t) can be directly measured by the technique of muon-spin rotation[27]. On the 
other hand, while the initial state preparation is not realizable, the interest lies in the 
equilibrium correlation function C (t) related to the cross-section of neutron scattering 
and the susceptibility x(o;)[T]. 

The effect of a harmonic environment is characterized by a spectral density 
J (no) = 2auj 1 ~ s uj s 6{uj c — u>) with the dimensionless coupling strength a, the upper 
cutoff uj c and the step function 6{x). An additional energy scale uj s is introduced, but 
only the combination au l s ~ s has fundamental significance. In this paper we assume 
that the high energy cutoff uj c is much larger than all other scales. It is convenient 
to set u s = u c /100 in this paper except the special notation. The index s accounts 
for various physical situations around the TSS. For example, in solid materials where 
acoustic phonon provides the most efficient damping mechanism, according to the 
Debye model, s = 3, which belongs to a super-Ohmic bath(s > 1). s = 1 and 
s < 1 stand for Ohmic bath and sub-Ohmic one, respectively. There are suggestions 
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to model 1/f noise found in experiments by a limiting sub-Ohmic bath|28] [29j [30] . 
In terms of the renormalization group, sub-Ohmic coupling represents a relevant 
perturbation [3 lj. Moreover, there are various claims that the particle is always lo- 
calized in the sub-Ohmic case for zero temperature based on noninteracting blip 
approximation (NIBA)[U [2]. As was pointed out in Ref. [6], the NIBA is unreli- 
able for studying the long-time behaviour of equilibrium correlation functions since 
it gives an exponential decay instead of an algebraic decay. Another tool for treating 
spin-boson problem, adiabatic renormalization, is invalid in the sub-Ohmic case[6]. 
Due to technical difficulties, little result about long time dynamics of s < 1 at zero 
temperature is known. 

The study of SBM with a sub-Ohmic bath for zero bias case has attracted much 
attention recently [5-12]. The flow equation method has been applied by S. Kehrein, 
A. Mielke and T. Stauber to this model [HI [TJ, [8] . Some important properties were pre- 
dicted, such as a transition from derealization for weak coupling to localization for 
strong coupling and the calculation on the equilibrium spectral function Ciui). How- 
ever, as was pointed out by Stauber, this approach did not yield a correct normaliza- 
tion condition nor a satisfying result about the Shiba relation for a certain parameter 
regime[8]. This approach seems difficult to calculate the quantum dynamics of SBM, 
especially the non-equilibrium correlation function P(t). Besides, the NRG method, 
a powerful numerical tool, was employed by Bulla, Tong, Vojta, et.al to investigate 
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the sub-Ohmic case and provide some reliable results for both static and dynamic 
quantities in the whole range of model parameters and temperatures [9], [Hll El EE]. 
They focused on the quantum critical behavior from a delocalized state to a localized 
one and found that the Shiba relation was fulfilled within an error of about 10% 
in the Ohmic case, but did not provide a check in the sub-Ohmic case[9]. Very re- 
cently, the NRG study for sub-Ohmic case shows that the dynamical properties of 
the delocalized phase are not dominated only by a single energy scale A r [12]. In their 
treatment, the energy scale ui s was set equal to the high energy cutoff u c . To the best 
of our knowledge, the dynamics transition between coherent oscillation and incoher- 
ent relaxation has not been given explicitly so far. Recently, Chin and Turlakov have 
used the variational method originally proposed by Silbey and Harris for the sub- 
Ohmic bath|5j. They also paid more attention to the transition between delocalized 
phase(the effective tunneling A r > 0) and localized one ( A r = 0) at both T = and 
T > 0. Even if they pointed out that dynamical and thermodynamics criteria for 
the transition should be expected different and sensitive to non-adiabatic mode. In 
their discussion, they have ignored the dynamical effect of the perturbation. Thus, 
it is difficult for the variational ansatz to make detailed statement about quantum 
dynamics beyond Born-Markov approximation. Moreover, they worked with u s ^ u c 
and io c — > oo. As far as we know, there is no numerical and analytical approach 
that can give the dynamical transition of nonequilibrium correlation function, exact 

5 



Shiba's relation and sum rule in the sub-Ohmic case. In the present paper, we ex- 
tend a unitary transformation proposed by one of our authors to investigate quantum 
dynamics of the sub-Ohmic SBM at T = 0[14j. 

The paper is organized as follows. In Sec. II, the hamiltonian is separated into 
the unperturbed part and the perturbed one based on a unitary transformation and 
perturbation theory. In Sec. Ill A, the localization-delocalization transition point will 
be clarified, compared with the results of flow equations and NRG. The calculations 
of quantum dynamics P(t) and C(t) will be given explicitly in Sec. Ill B. Then, It is 
verified that real and imaginary parts of x{ u ) satisfy the Shiba's relation analytically 
and numerically in Sec. Ill C. The coherent-incoherent transition is also discussed in 
detail in this part. Finally, Sec. IV gives discussion and further analysis of our results 
and approach . Usually, people believe that perturbation approach is not good for the 
dissipative SBM because of the infrared divergence in calculating the renormalized 
tunneling frequency and other physical quantities by perturbation expansion. Here 
we attempt to get rid of the divergence by using a unitary transformation. When 
one take the scaling limit A << uo c from our obtained results, it will reproduce the 
known or exact one. This approach works well for the low-temperature region and 
weak coupling case with < A < uj c far away from the scaling limit. 
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II. Unitary transformation 



In order to take into account the correlation between spin and bosons, we present 
a treatment based on the unitary transformation to H: H' = e s He~ s , where the 
generator of the transformation is 

S = J2^M-h)a z . (2) 
k ZuJk 

A A;— dependent function introduced in the transformation corresponds to the 
displacement of each boson mode due to the coupling to the TSS[H1 E2]- Its form 
will be determined later by the perturbation theory. 

The transformation can be done to the end and the result is H' = H' Q + H[ + H' 2 , 
where 

H' Q = - X - V Aa x + J2u k blb k - £ J-^(2 - (3) 

k k k 

V = eM-Y,^-M (4) 
k ZLJ k 

H[ = lj2 9k(l - tk)(A + h)a z - \ V Aia y £ - &»), (5) 

1 k 1 k U k 

H' 2 = ~Aa x (cosh{£ ^ k (b\ - b k )} - v ) 
~Aia y (sinh{£ ^Ub\ ~ h)} -^—^A ~ 

1 \ k Uk k U k J 



(6) 



H' is the unperturbed part of H'. Obviously, since the spin and bosons are decoupled 
in this part, H' can be solved exactly. The eigenstate of H' is a direct product: 



|s)|{rifc}), where \s) is the eigenstate of a x : \s\) = y ^ j or |s 2 ) — y ^ j, and 
| {rife}) is the eigenstate of bosons with n k bosons for mode k. In particular, |{0fc}) is 
the vacuum state in which n k = for every k. The ground state of H' is 

l»b> = l«i>|{0*}>. (7) 

//{(first-order terms) and H' 2 (including second order terms and higher ones )are 
treated as perturbation and they should be as small as possible. For this purpose 
is determined as 



& = ^-r- (8) 
Substituting this form into Eq. (5), one has 

H 'i = h A E — r-A [b ^ az ~ l<yy) + bk{(Jz + l(Jy)] - (9) 

2 k uj k + f]A 

It is easy to check that H[\g ) = 0. Thus, by choosing the form of the matrix 
elements between the ground and the lowest excited states are zero (we show them 
in the following). It will become possible to take a perturbation treatment based on 
the division of the transformed hamiltonian. Besides, in transformed hamiltonian H' 
7]A is the probability of diagonal transition of bosons which describes the coherent 
tunnelling motion of particle, rj is determined in Eq.(4) to make Tr(p B H' 2 ) = 0, 



where ps = exp(— /3ifs)/Trexp(— (3Hb) is the equilibrium density operator of bosons 
(Hb = J^k^kbl-bk)- Thus one can see that the tunneling is renormalized by this factor 
r) arising due to the dressing of bosons coupled with the TSS. In other words, the 
particle is surrounded by bosons cloud as it tunnels between the two states. 

It is noticeable that < < 1, which determines the intensity of the correla- 
tion between spin or particle presentive for subsystem and bosons in bath: ~ 1 
if the boson frequency (uik) is much larger than the renormalized tunneling of sub- 
system (rjA) while ((. C 1 for Uk <C ?]A. Since the transformation in Eq.(T5]) is 
essentially a displacement one, physically, one can see that high-frequency bosons 
(ujk > T]A) follow the tunneling particle adiabatically(instantaneously) because the 
displacement is g^k/^k ~ gk/^k, which leads to a dressed particle. On the other 
hand, low-frequency bosons uJk < f]A, in general, are not always in equilibrium with 
the tunneling particle, and hence the particle moves in a retarded potential arising 
from the low- frequency modes. When the nonadiabatic effect dominates, ojk <§C rjA, 
the displacement g^k/^k ~ gk/v^ <C 1, is substantially reduced. Therefore, the 
effects of each boson mode on subsystem are in nature treated separately by 

The lowest excited states are |s2)|{0fc}) and |si)|l&), where is the number 
state with = 1 but ny = for all k! ^ k. It is easy to check that {go\H' 2 \g ) = 
(because of the form of r\ in Eq. (4)), ({OfcH^I-f^lfi'o) = 0, (1^1 (siji^ls'o) = 0, and 
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({0 k }\(s2\H' 2 \ Sl )\l k ) = 0. Moreover, since H[\g ) = 0, we have ({0 k }\ (s 2 \H[\g ) = 
and (lfc|(si|iif{|<7o) = 0. Thus, we can diagonalize the lowest excited states of H' as 

H' = — -r]A\g ) (g \ + ^ E\E)(E\ + terms with higher excited states. (10) 
2 E 

The diagonalization is achieved through the following transformation [24J: 



\s 2 )\{0 k }) = Y,x(E)\E), (11) 

E 

\ Sl )\l k ) = Y,Vk{E)\E), (12) 

E 

\E) = x(E)\s 2 )\{0 k }) +J2yk{E)\ Sl )\l k ), (13) 

k 



where 



, ~ 1/2 

r[t) ' ' ^ (E + V A/2 - cu k )\ 



V 2 



(14) 



y^ E ) = f _l a/9 — (15) 

with Vfc = r]Ag k £ k /ui k . E's are the diagonalized excitation energy and they are 
solutions of the equation 

E-^-V S = 0. (16) 

III. Quantum dynamics 
A. Delocaliztion 



The Ohmic spin-boson model has nontrivial dynamics only for a < 1 when A <C c<j c . 
Scaling arguments, flow equations and other methods give a renormalized tunneling, 
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A r = A(^-) 1_r ^[Tl [15] . When a goes to 1, one gets A r = which is a transition 
from derealization for weak coupling to localization for strong coupling. However, 
for the sub-Ohmic bath s < 1, there is some confusion about the existence of such 
a transition within some various approximation approaches. On the one hand, the 
TSS is localized for nonzero sub-Ohmic coupling at T = which is predicted by 
the NIBAJTJ [2J. On the other hand, from the mapping of SBM to Ising model, the 
corresponding results turn out a transition as a function of coupling [6], [15]. Moreover, 
the results of flow equations and numerical renormalization group also support it [HI [9]. 
However, the transition point is less clear for general A and s. 

The renormalization of tunneling can be calculated as 

-"fA ( " A ' r ' + ^--h-^^r 1 } ■ w 

Here uj' s = uj s /uj c and A' = A/u c . From this equation, it is self-consistent to deter- 
mine rj. It is clear to see that the eigenstates of the unperturbed part in H' are of 
superposition of the TSS if r\ > 0, whilst TSS becomes localized if r\ = 0, and no 
tunneling is possible. According to the analogical treatment of the transition con- 
dition in the work of Kohrein and Mielke[6], the critical transition condition from 
delocaliztion to localization can be derived. 

/ , 

r] = exp < — aoo s 
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>\s-l 



TCS 



smns 



+ E 

n=0 



n+1 

1 

s — n — 1 



;is) 



For the scaling limit A' <C 1, one gets 

TVS 1 



1] = exp < — auj'g' 



sin its s — 1 



(19) 



If the condition 



auf- < e- 1 (20) 

s sm(7rs)A /1 - s _ v ; 

is satisfied, the solution for r\ is finite and satisfies 

otherwise = 0. It can be seen that, in the sub-Ohmic case, there exists a quantum 
transition boundary separating a delocalized phase for a < oti from a localized phase 
for a > oil. In the delocalized region, the renormalized tunneling between the two 
states is finite, whereas it is renormalized to zero in the localized phase(?7 = 0). When 
s — > 1, from the condition (!20l) . it is easy to get a\ = 1 in the scaling limit, which 
agrees with the known results. Our results provide a direct evidence for a phase 
transition for all < s < 1. Furthermore, the transition point ai is continuous as a 
function of s from the sub-Ohmic to Ohmic case, which is in good agreement with 
the NRG results [9l [TO]. Compared with the transition condition Eq. (16) in Ref.[6] 
obtained by flow equations, it seems that there is a discontinuous behavior when s 
goes below to 1. Moreover, no transition happens for s > 1, in other words, the 
system is always delocalized in the super-Ohmic case. As is shown in the condition 
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(|20|) . the critical coupling a\ follows a power law as a function of the bare tunneling, 
ai oc A 1_s , which is in good agreement with the conjecture of NRG [TO]. From above 
discussion, it indicates that the transition condition is sensitive to the index s in the 
sub-Ohmic case. 

f] as a function of the dimensionless coupling strength is shown in Fig. 1. 

There is a discontinuous jump from a finite value to rj = as a — > a\ for < s < 1, 
namely, it is a first-order transition which is different from the NRG result [HI [TO] . 
However, for s — 1, r\ can continuously change from one to zero as a increases from 
zero to oli ~ 1. For a > 1 in the scaling limit, the delocalized transition occurs which 
agrees with the known result in literature. 

Our approach has given a clear evidence for the transition from a delocalized 
phase to a localized one in the sub-Ohmic case. The phase boundaries are shown in 
Fig. 2(a) which are determined by the vanishing of renormalized tunneling. NRG 
results are shown for comparison. It is seen that our data agree well with those of 
NRG for s > 0.6. At the same time, there is some deviation from our data for s < 0.5 
due to the NRG discretization [TD] . As displayed in Fig. 2(b), the critical coupling 
o>i follows a power law as a function of the tunneling, a\ oc A 7l ~ s for A << uj c (see 
Eq. [20]), also predicted by the NRG methodpJJ]. NRG data are provided for a fair 
comparison in Fig. 2(b), too. The calculated transition points are in agreement 
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with those obtained by the NRG method. However, our result predicts a first-order 
transition, in contrast with the second-order transition predicted by the NRG. 

B. Correlation function 



The main interest of quantum dynamics is the non-equilibrium correlation P(t) and 
the symmetrized equilibrium correlation C(t). When the initial state can be prepared, 
the evolution of the state is of primary interest which can be described by P(t). 
P(t) = (b,+l\(+l\e iHt a z e- iHt \ + 1)|6,+1) is defined in the Ref. [T], where | + 1) is 
the state of bosons adjusted to the state of a z = +1. 

Because of the unitary transformation (e s a z e~ s = a z ) 

P(t) = ({O k }\(+l\e iH,t a z e- iH,t \ + l)|{0 fc }>, (22) 

since e | + l)\b, +1) = | + l)|{0fc}). Using equations(10-16) the result is 

P(t) = - Y x 2 (E) exp[-i(E + rjA/2)t] + - Y X \ E ) exp[*(£ + r]A/2)t] 

^ E ^ E 

= — f due- 4ut (u-riA-y i ^ ) 

4m Jc \ fc ^ + iQ -Uk) 



+— / due™ (u-TtA-Y S ) , (23) 



47ri Jc \ ^ uj — i0+ — uJk / 

where a change of the variable uj = E + 77 A/2 is made. The real and imaginary parts 
of J2k V£ I i u ^ ~~ u k) are denoted as R(uj) and =py(a;), respectively. 

R{u) = V(^^-) 2 ^- 

V ^k + vA J uJ-uj k 
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-2a^-°u c ( V Ar *** (24) 

Jo (x — Lu')[x + riA'Y 



(x - LO')(x + TjA'y 

7 (w) = 2«vr^ 1 -U( Z7 ^ 7 ) 2 (^) s 0(^ - «), (25) 

where u' = uj/uj c for simplicity. For general s < 1 the integration in Eq. ([2lj) should 
be done by Residue Theorem, 

;1 _ s 2 f " (1 - g )(-^AQ"+ 2 + u's(-7]A') n+1 - (o/)"+ 2 
Mom = — 2aoo„ ujA — ) < > 



(r/A' + uj')r)A' 7r 



(s - l)(r/A') s + s(t/A') s ~V + (u/) s cos(tts) } .(26) 



1 + rjA' sin(7rs) 

The integral in Eq. (l23il can proceed by calculating the residue of integrand and the 
result is P(t) = cos(u;ot) exp(— jt), where u is the solution of equation 

u - r]A - R(u) = (27) 

and 7 = 7O7A) = aTTul^ s (T]A) s /2 (Wigner-Weisskopf approximation) [33] . The be- 
havior of P(t) is of the form of damped oscillation. One can prove that the solution 
uj of Eq. (1271) is real for small coupling. As the coupling a increases, the solution 
ujq becomes imaginary, thus P(t) demonstrates the incoherent dynamics. As a result, 
there exists a critical point corresponding to the coherent-incoherent transition. In 
other words, for the critical case, one can have ujq = and P(t) = exp(— 7 c t). For 
a > a c , one has P(t) > for all times. Meanwhile, the behavior of damped oscilla- 
tions disappears and that of pure incoherence displays. When A' 1, from Eq. [2] 
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one gets 



sin(7rs) ,T]A' 



l-s 



(28) 



2vr(l- S ) v cu> 



For s = 1, one gets a c 



2 which is the same as known before. It is clear to see that 



a c is also proportional to A 1 s for small A. It turns out the sensitivity of the critical 
coupling to bath structure. 

Fig. 3(a) shows the time evolution of the non-equilibrium correlation P(t) for the 
fixed coupling au'}~ s = 0.1 and A/u s = 10 with different bath types. The damped 
oscillation exists when a < a c , while it decays fast for a ~ a c . That is to say, 
the tunneling regime a < ai, namely rj > 0, consists of two qualitatively different 
regions, distinguished by the presence (for < a < a c ) or absence (for a c < a < a>i) 
of tunneling oscillations in quantum dynamical quantities j2TJ. It is predicted that 
there is a boundary between the phase coherent region and pure incoherent one( the 
exponential decay), which is shown in Fig. 7(a). 

The behavior of P(t) is shown in Fig. 3(b) for s = 0.9 and a = 0.1 with different 
tunneling A/u s = 1,5,10,20, and 30, respectively. The less A is, the faster P(t) 
decays. It indicates that there is a nonscaling behavior for dynamical quantities in 
the sub-Ohmic bath in contrast with a scaling one in the ohmic case. 

Since e s a z e~ s = a z , the retarded Green's function is 



G{t) 



i9(t) {[exp(iH't)a z exp(-iH't), a z }+)' , 



(29) 
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where (...)' means the average with thermodynamic probability exp(—/3H'). The 
Fourier transformation of G(t) is denoted as G(u), which satisfies an infinite chain 
of equation of mot ion [31] . We have made the cutoff approximation for the equation 
chain at the second order of g k and the solution at T = is 



G(U) ui-vA-Zk V k 2 /(u - cu k ) + u + V A - Ek V k 2 /(iO + cu k ) ■ m 



The equilibrium correlation function 



C(t) = ^Ti{exp(-(3H)[a z (t)a z + a z a z (t)}}/TT[exp(~(3H)] 



1 r°° j3u 

— I <iu;coth( — )lm.G{u) exp(— iut) 
In J-oo 2 



2tt 

1 f°° li^) 

= — duj- — -— -cosiut). (31) 

n Jo [u - r]A - R(uj)} 2 + 7(cj) 2 v 1 y 1 

For general value of a < a c , C(t) may contain both terms of the exponential decay 
ones and the algebraic decay ones. For the special value a = a c , the exponential 
decay terms disappear. In the long-time limit the first non-zero algebraic decay term 
dominates which is ~ — l/t s+1 . The time evolution of C(t) is shown in Fig. 4 for 
s = 0.8, A/u s = 0.1 with different couplings a = 0.1,0.3 and 0.5, respectively. 
The time evolution of P(t) is also provided to a fair comparison in this figure. It is 
remarkable that there are distinguishable characters between C(t) and P(t) for larger 
coupling, while the difference between them is not apparent for small coupling. This 
is an indication that P(t) and C(t) possess a similar structure for small coupling at 
intermediate times. However, they differ distinctively at long time limit. 
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In Fig. 5 we show C{uo) for s = 0.6 (Fig. 5a)and s = 0.3(Fig. 5b). For 
weak coupling, there is a peak near the renormalized tunnelling A r . However, with 
increasing coupling, it is observed that part of the spectral weight has transferred to 
lower frequencies with a shoulder feature. It indicates that the dynamical properties 
in the sub-Ohmic case is not determined only by a single energy scale A r which agree 
with the conclusion of the NRG [32] 

C. Shiba's relation and sum rule 

The susceptibility xfa) = ~~ and its imaginary part is 

X y ' [u-rjA- R{uj)Y + 7 2 (u;) [uj + 77A + R{-u;)} 2 + 7 2 (-cu) ' 1 J 
The uj -f limit of S(cu) = x"(u)/J(u) is 

X"{U) 7T 

JM ~ (r / A + J R(^ = 0)) 2 - (33) 

With a Kramers- Kronig relation and a fluctuation-dissipation theorem, the static 
susceptibility Xo can be directly extracted 

* = »)=if^. (34) 
2 jo 

Thus, there exist an important relation connecting the zero frequency behavior of the 
spectral function to static susceptibility, which is known as the Shiba's relation [T8| 
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UHl [201 [231 EE 



lim 



C{u) 



1. 



(35) 



It constitutes an important check of our approach. 

Approximation schemes such as NIBA and numerical approaches based on Monte 
Carlo simulations can not be used to verify the Shiba's relation since they failed to 
predict the long time behavior [7]. Results from both real time renormalization group 
and NRG in the Ohmic case have a maximal relative error of around 10%[9l[T9]. Flow 
equation approach based on infinitesimal transformations gives its result for the sub- 
Ohmic case with error 10% or so [8]. In Table I, the generalized Shiba relation is 
verified for various couplings, tunneling matrix and bath types. It shows that the 
Shiba's relation is exactly satisfied for a < a c in both the sub-Ohmic and Ohmic 
cases. Outside the range, the agreement is still good but no longer exact. In the 
present paper, we consider carefully the behavior of coherent dynamics for a < a c . 
For s — 1, the Shiba's relation has been proved analytically by our approach in Ref. 



The sum rule C(t = 0) = 1 is an important relation in the equilibrium correlation 
function which constitutes another significant check of our approach. The numerical 
solution of our approach can yield a satisfactory result, which is also shown obviously 
in Table I. As is told in the paper of StauberjS], the sum rule is not fulfilled by flow 



[14j. 
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equation method. Independent of the bath type and the cutoff oo c the sum rule is 
performed well in the present approach, while in the flow equations it only yield 80% 
or so. 

Since different dynamical quantities may be associated with different initial prepa- 
rations of the system, quantum coherence may be more or less sensitive to dissipation [25] 
Therefore, there are disparate critical values of the damping strength for different co- 
herence criterion. Considering the character of the spectral function S(u), one gives 
its appropriate coherence criterion. For weak damping, a < a*, the function S(u) 
exhibits two inelastic peaks at finite frequency u p ~ 77A. It can be seen that the 
width of peak 7(0;) is less than the value of uo p . At the critical damping a*, one has 
7(0;) = Up, and the two peaks merge into a single quasielastic peak centered at uj = 0. 
For instance, in the scaling limit, for s = 1 , the critical damping value is determined 

by 

r/A 7r^ ~^ i]A i]A rjA ' ^ 

and 

2«:(t-V) 2 = 1- (37) 

One gets a* = 0.325, which agrees well with a* = | or a* ~ 0.3 obtained by the 
renormalization group numerical results [9]and other conjectures [201 E2J- One can 
observe that the " quasiparticle peak" disappears at a = a* in Fig. 6 (b). Physically, 
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it indicates that the crossover between the underdamping oscillating behavior and 
overdamping decaying one occurs at a = 0.325 at zero temperature. 

In Fig. 6 the renormalized spectral function S(u)A 2 is shown for s = 0.8 (a) and 
s = 1 (b) with A/u s = 1 , u) s /oj c = 0.1 and various coupling strength a. S(u) has 
a double-peak structure for a < a*. Only the u > part is shown in Figs. 6(a) 
and 6(b). For a > a* c there is only one peak at uo = 0. The critical coupling a* 
corresponds to the crossover from underdamping to overdamping oscillations. We 
obtain a* = 0.34 for a subohmic bath with s = 0.8 which is in good agreement 
with a* ~ 0.3 obtained by the flow equations[7j. Fig. 6(c) shows the S(u) versus u 
relations for fixed a = 0.3 and s = 0.9 with A/cu s = 1, 5, 10 and 30, respectively. 

The dynamical phase diagrams are given in Fig. 7(a) (s = 1) and Fig. 7(b) 
(s = 0.8) for uj s = td c . In both figures, the transition boundary between the in- 
coherent relaxation and coherent oscillation is shown in a solid line, while that of 
underdamping-overdamping transition is shown in a dashed line. In the scaling limit 
(A' << 1) for s — 1, a c — 0.5 and a* = 0.325 are consistent with the known 
resultsjTJ[20l[22]. Further, in the Ohmic case, the coherent-incoherent transition cou- 
pling a c = |[1 + r]A/u c ]. With increasing tunneling, both a c and a* become larger, 
which presents a contrast with the tendency of the boundary obtained by Monte 
Carlo simulations |18j. Moreover, note that the critical coupling is sensitive to the 
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bath type, especially in the scaling limit. As displayed in Fig. 7(b), the critical cou- 
plings a* and a c nicely follow power laws as functions of the bare tunneling, a oc A 1_s 
for small A (in the case s = 0.8). Our fit data for A' < 0.01 are well consistent with 
the conclusion. It is pointed out by our approach that the critical couplings in ad- 
dition to the delocalized-localized transition one are always proportional to A 1_ ~ s for 
small A, which hints that the critical value is sensitively dependent on the structure 
of bath. 

IV. Summary and discussion 

The dynamics of SBM with a sub-Ohmic bath is studied by means of the perturbation 
approach based on a unitary transformation. Our approach is quite simple, whereas it 
correctly gives P(t), C(t), Shiba's relation, and reproduces nearly all results of previ- 
ous authors that we note. By means of our approach, the dynamical transition point 
a c is calculated for the first time. The main results include: (1) the non-equilibrium 
correlation P(t), the susceptibility x"( w ) an d the equilibrium correlation C(t) are 
analytically obtained for the general finite A/u c case; (2) a critical transition point 
(ai) from derealization to localization exists in the sub-Ohmic case; (3) for a fixed 
tunnelling, as the coupling with environment increases, a crossover from coherent to 
incoherent tunnelling appears. In other words, for a > a c , the coherent dynamics 
disappear; (4) there is a critical transition point a* from underdamping oscillation 
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to overdamping one in the susceptibility When a < a*, the spectrum of 

susceptibility is of a double-peak structure, while there is only one peak at u> = for 
a > a*; (5) the Shiba's relation and sum rules are exactly satisfied for a < a c ; (6) 
the dynamical phase diagrams in both the sub-Ohmic and Ohmic cases are mapped 
out. For small tunneling, all of three critical coupling a c , a* are proportional to 
A 1 - 3 . 

Our treatment is essentially a perturbation one. The transformed hamiltonian is 
divided into the unperturbed part H' and perturbation H[ and H' 2 . The unperturbed 
part gives the intrinsic coherent tunnelling motion with diagonal transition of bosons. 
H[ is related to the single-boson non-diagonal transition and H' 2 to the multi-boson 
transition. In our perturbation treatment we take into account only the single-boson 
non-diagonal transition (H[) which plays an important role for weak coupling and 
lower temperature. H' 2 with multi-boson non-diagonal transition is neglected in our 
treatment. At strong coupling and high temperature, the multi-boson non-diagonal 
transition may play an important role. Since our perturbation treatment keeps the 
contribution of the single-boson non-diagonal transition and drops that of the multi- 
boson non-diagonal transition, we can not make definite statement about the exact 
nature of the localization-delocalization transition. Our perturbation treatment could 
not be applied to study the localized phase beyond the critical point of the sub- 
Ohmic case(o; > on). Thus, the approach can not predict the oscillation behaviors 
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on short times in the localized phase for strongly sub-Ohmic case. The recent NRG's 
results have captured the coherent oscillation even in the localized phase for small s 
and demonstrated it explicitly [12J. In principle, our approximation can be applied 
satisfactorily for a < a c , whereas the approximation might fail when a ~ ai. 

At lower temperature and weak coupling, the diagonal transition of bosons dom- 
inates with coherent tunnelling motion. As the coupling increases, single-boson and 
multi-boson non-diagonal transitions come to prevail against the diagonal transition, 
which result in decoherence. Our perturbation treatment keeps the contribution of the 
single-boson non-diagonal transition and drops that of the multi-boson non-diagonal 
transition. Based on the perturbation treatment, our approach can reproduce the 
well-known results for the Ohmic case, such as ai = 1 and a c = 1/2. Besides, Our 
result also turns out that the dynamical property is not determined only by the single 
energy scale r\ A. On the one hand, a non-scaling behavior of dynamical quantities 
is shown in Fig. 3(b). On the other hand, from the results of C(u) for s = 0.6 and 
s = 0.3 (Figs. 5(a) and 5(b)), one can see that a pronounced shoulder in C(u>) comes 
to appear with increasing coupling which agrees with the NRG results. These com- 
plicated behaviors result from the non-diagonal transition, which is the contribution 
of H[ in our perturbation treatment. At the same time, Shiba's relation and sum rule 
are exactly satisfied for a < a c in both the Ohmic case and sub-Ohmic one ( even for 
s<l )which constitute important checks of the approach. In general, this method 

24 



works well for weak coupling with a < a c . Thus, this approach can be applied in the 
sub-Ohmic case to study the coherent-incoherent transition. 

The function ^ plays an important role in our approach. It eliminates the in- 
frared catastrophe completely. In conventional perturbation theory, from the original 
Hamiltonian H, the dimensionless expansion parameter is g\joj\. Thus, the renor- 
malized tunneling becomes zero, and the system is always localized because 



For Ohmic bath s = 1 the integrand becomes 2a /u which is logarithmic divergent in 
the infrared limit. There also exists a low-energy divergence (^rj) for a sub-Ohmic 
bath s < 1. On the contrary, for the coupling in transformed Hamiltonian, H[, the 
expansion parameter is gl^l/ojl ~ 2auj s uj l s ~ s / (u + rjA) 2 , which is finite in the infrared 
limit for s > 0. Therefore, from the viewpoint of mathematics, the introduction of 
gets rid of the divergence. Physically, the disappearance of divergence arises due 
to the energy scale separation by with a special treatment of the low frequency, 
nonadiabatic modes. 

The low frequency behavior of the spectrum density function in SBM model de- 
termines the long time behavior of TSS. All quantum dynamics properties are very 
sensitive to the low energy part in spectral structure, especially in the sub-Ohmic 
case. On the other hand, in the high frequency limit (cj^ >> A), bath modes follow 




(38) 
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instantaneously to the tunneling motions, namely, the displacement of boson due to 
the coupling to subsystem is large, whereas near the low frequency limit (u>k << A), 
nonadiabatic modes couple weakly to the subsystem with coupling strength gk^%, its 
displacement is small. Thus, while all boson modes are treated by the function 
their different contribution to the dressed TSS has been distinguished with respect 
to the scale of boson energy. Therefore, even if our approach is, in principle, a per- 
turbation one, those reasons above give the intrinsic base of our present results, such 
as P(t), C(t), the dynamical transition points a c and a* c and exact Shiba's relation. 
When s — > 1, our results agree well with exact or known ones. Besides, our approach 
can be straightforwardly extended to other more complicated coupling systems. 
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Figure Captions 



Fig.l rj as a function of the dimensionless coupling strength a(^) 1 s for different 
bath types. 

Fig. 2 (a)The delocalized-localized transition point ai as functions of s. (b) A 
dependence of the critical coupling a\ for various bath types s. The NRG data are 
also shown for comparison. NRG parameters here are A = 2, Nb = 8, and Ns = 100. 

Fig. 3 (a)The time evolution (A r t = r]At)oi the non-equilibrium correlation P(t) 
for auj's~ s = 0.1 and A/uj s = 10 with different bath types s. (b) P(t) for s = 0.9 and 
a = 0.2 with various values of A. The inset shows P(t) for a = 0.1 with different A 
in the Ohmic case. 

Fig. 4 The time evolution of P{t) and C{t) for s = 0.8 and A/uj s = 5 with 
different couplingsa = 0.1,0.3, and 0.5, respectively. Note that a c = 0.53718. 

Fig. 5 The equilibrium correlation function C(uo) for (a) s = 0.6, A' = 0.1 and 
au)'}~ s = 0.05, 0.08, 0.12, and 0.13; (b) s = 0.3 with various values of a, and A' = 0.1 
(w s = u c = 1). we find a c = 0.1271 for s = 0.6 and a c = 0.0327 for s = 0.3. 

Fig. 6 The renormalized spectral function S(uj)A 2 as a function of oo for s = 0.8 
(a) and s = 1 (b) with A/uo s = 1 , uo s /uo c = 0.1 and different a. (c) S(uo)A 2 as a 
function of uo for s = 0.9, a = 0.3 with A/uo s = 1, 5, 10, and 30, respectively. 

Fig. 7 The dynamical phase diagrams for s = 1(a) and s = 0.8(b) with oo s = u c . 
The solid lines are the transition boundaries between the coherent oscillation and in- 
coherent decay. The dashed lines stand for the boundaries between the underdamping 
oscillation and overdamping one in coherent region. The dotted lines in fig6. (b) are 
fits to the a oc A 1_s law using the A < 10~ 2 points only. 
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Table Captions 



TABLE I: Representative results from the numerical solution with parameters 
chosen by the spectral density J(ou) = 2aujl~ s u) s O(ui c — uo) and with the controlling 
precision 10~ 5 for iteration. R = [lim^o C(^)/^(^)]/(2xo) 2 -The numeric error for 
the Shiba relation and sum rule is at least less than 10~ 6 and can be improved by 
increasing the accuracy of numerical calculations. 
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